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ABSTRACT 

A free-floating wave meter was used to obtain a time 
series of wave-induced pressure variations in deep water. 

The meter’s response was roughly checked against a piston- 
type pressure testing apparatus. 

Techniques in the operation of the meter were developed 
involving a slight modification in the suspension of the 
sensor. 

Power spectra of the pressure variations were obtained 
from the digitized wave record via the amplitude spectra 
given by the fast Fourier transform. 

Characteristics of the power spectra of the metered 
waves were used to evaluate the meter by comparing them with 
the waves that might be anticipated (on the basis of 
theory) from the existing surface wind field and with visual 
observations of the metered waves. 
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I. OBJECTIVES 


The U.S, Naval Postgraduate School acquired a free 
floating wave meter designed to measure deep water wind 
waves in the sea over the band of frequencies associated 
with swell and locally wind generated waves. It is desire- 
able to make operational use of an instrument of this kind. 
Most wave metering requires a physically stable reference 
platform from which to record deep water wave activity, 
or it requires that the reference platform's motion be 
known. These tedious and expensive requirements are avoided 
with the free floating meter. In order to facilitate the 
use of this instrument in future research, the objectives 
of this paper are: 

1. To estimate the response of the instrument by a 
direct calibration procedure; 

2. To providB a detailed account of techniques and 
equipment used in making wave measurements at sea with 
the instrument; 

3. To analyze the wave records and to compare the 
results with independent estimates of the wave parameters 
made for the same wave system. 


9 



















































































II. THEORY 


The free floating wave meter has a pressure sensor hung 
below a float which moves vertically with the free water 
surface. 

The wave pressure for wave number, k, of the fixed depth 
is given by: 

-kz^ 

p^(k; z^,t) = p^(k;t)e 

where p^ is the static (or slowly time varying) component of 
the pressure due to the rise and fall of sea level. 

For the bobbing meter suspended at the nominal (mean) 
depth, z^, below the undisturbed sea surface, there are two 
components to the pressure variation: the dominant contri¬ 
bution is that due to the rise and fall of the sensor as the 
surface float is displaced, this being Identically -p^(k,t), 
and the other is the change in environmental pressure p^. 
Thus, the total pressure change experienced by the meter is 
approximately: 

-kz-. 

p^(k;z^,t) = -p^(k;t) + p^(k;t)e 
-kz. 

p^(k;z^,t) = -p^(l-e 

Therefore, assuming that the sensor hangs vertically from 
the float, the instrument will detect pressure fluctuations 
corresponding to the vertical motion of the float. 

If one wants to measure the sea surface displacement as 
indicated by the bobbing motion, the meter must be suspended 
at a sufficient depth z^, such that: 
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p^(k;z^,t)_ = -p^(k;t) 

That is, accurately describes the surface wave action if 

-kz -k z 

e -'-«1 for wave numbers k^k^. Letting e o = .1 and 

z = 300 feet (the depth of suspension for the sensor in the 

measurement made in this study), pressure measurements from 

waves with wave numbers greater than .0077 feet"^, or 

frequencies greater than .079 seconds”^ will not be in error 

by more than 10^. Gaul and Brown [1] give more detail in 

the theory of a free floating wave buoy. 
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III. DESCRIPTION OF THE INSTRUMENT SYSTEM 


The instrument system has a transducer-buoy package, 
cc signal conditioning deck unit, and recording units 
CFlEure 1). 

The transducer-buoy package (Figure 2) consists of a 
sensitive pressure transducer, small diameter (.060 in.) 
sea cable, and buoy hull, which contains a modulator, 
r^dio transmitter, batteries, and antenna. The transducer 
detects pressure changes by means of a vibrating wire whose 
fiiequency varies directly with the applied hydrostatic 
pressure. The sensitivity and frequency response are 
discussed in Section IV. The transducer and buoy electronics 
convert frequency changes into an FM signal with a center 
frequency of approximately 27 MHz. The components of this 
package are designed to minimize water and wind drag. The 
transducer has a streamlined torpedo shape of length 30 cm. 
and diameter 3 cm. The buoy is circular with a 60 cm. 
diameter in horizontal section, and has an elllptlcally 
shaped vertical cross-section of depth 1^ cm. 

The signal conditioning deck unit consists of a modified 
citizen’s band radio receiver, a discriminator, a d.c. 
amplifier, and a band pass amplifier. The deck unit condi¬ 
tions the signal from the buoy so that a ten foot peak-to- 
peak wave results in a five volt peak-to-peak signal. The 
output signal does not represent absolute pressures; it 
ref-lects changes from the static pressure. 
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Figure 1. Instrument System 




































































































































Figure 2. Transducer-buoy Package 




































































The recording unit provided with the instrument is a 
Leeds & Northrup Speedomax H chart recorder. In addition 
to the chart recorder, the output signal from the deck 
unit was sent through a d.c. amplifier and recorded on a 
Sangamo 3500 l^-track magnetic tape recorder. A signal 
generator was also used to put a 60 Hertz reference signal 
on a track other than the data track during the period of 
recording in order to identify the data to be processed. 

A double-beamed oscilloscope was linked to the playback 
outlets of the tape recorder in order to verify that the 
data and reference signal were being recorded. 

More detailed Information about the instrument can be 
found in Gaul and Brown [2]. 
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IV. CHARACTERISTICS OF THE SYSTEM 


This section discusses four characteristics of the 
instrument: frequency measurement limitations due to size 

of float, calibration dial characteristics, frequency 
response of the transducer-buoy electronics, and the 
response factor of the instrument, 

A. FREQUENCY MEASUREMENT LIMITATIONS DUE TO SIZE OF FLOAT 
The accuracy of pressure measurements, in part, depends 
upon the ability of the buoy to follow the surface of a wave. 
The buoy does not respond to waves with a length less than 
twice the diameter of the float (2,33 feet). Therefore, let 
the minimum observable v/ave length, L^, equal ^.66 feet, 
and utilize the deep water equation 



where T is the period of the wave in seconds and L^ is in 
feet. From this equation, T is 0.95^ seconds or frequency 


is 1.05 sec”^. Therefore, theoretically, the buoy does 


not respond fully to passing waves with frequencies greater 


than 1.05 sec ^ (i.e. ^ I sec. waves). 


B. CALIBRATION DIAL CHARACTERISTICS 

The deck unit has a calibration dial used to maintain 
a 5 mvolt peak-to-peak output signal, which is equivalent 
to full scale deflection on the chart recorder for the same 
maximum wave height v/hen the transducer is used at various 
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depths. For example, a dial setting of 500 when the trans¬ 
ducer is at a depth of 300 feet gives a 5 mvolt peak-to-peak 
signal when a ten foot peak-to-peak wave passes. Figure 3 
shows the dial setting for various depths to get full scale 
deflection for 10 foot waves. Figure ^ and 5 show dial 
settings vs. peak-to-peak voltages at constant depths for 
10 foot waves. 

C. FREQUENCY RESPONSE OF THE TRANSDUCER-BUOY ELECTRONICS 

Table I on frequency response of the transducer-buoy 

electronics was obtained from the manufacturer of the system 
Bissett-Berman Corporation, Hytech Marine Products. The 
table Indicates that the pressure transducer is more respon¬ 
sive to pressure changes as the absolute pressure Increases. 
By applying Increments of pressure at various absolute 
pressures with a piston-type pressure gauge tester, it was 
concluded that the Instrument does not record the magnitude 
of the pressure changes with sufficient sensitivity unless 
the absolute pressure is greater than ^0 to 60 psi. There¬ 
fore, the Instrument does not function properly in water 
with a depth less than 100 feet. This characteristic made 
it Impossible to calibrate the meter by direct comparison 
with meters placed at near shore locations. 

D. RESPONSE FACTOR OF THE INSTRUMENT 

An attempt to compute the response factor of this Instru 
ment was made by obtaining power spectra of a known input 
and its response. The piston-type pressure gauge tester was 
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Figure 3i Depth vs. Dial Setting for 10 foot Wave 
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Figure 4. Output Voltage vs. Dial Setting at 300 feet for a 10 foot 
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Figure 5. Output Voltage vs. Dial Setting at 150 feet for a 10 foot Wave 
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FREQUENCY RESPONSE OF TRANSDUCER-BUOY ELECTRONICS 





















used to apply a 2 psi Impulse over 2 seconds. After the 
impulse, the pressure was held steady for approximately 20 
seconds. 

Three separate trials were made for which the input, as 
constructed from visual readings of the pressure gauge and 
stopwatch dials, is essentially as shown in Figure 6. The 
power spectrum analysis of this input and the analyses from 
the three outputs are compared in Figure 7. Variations in 
the spectra of these three responses were due to the dif¬ 
ficulties involved in creating three identical pressure 
impulses with the pressure gauge tester. 


Pressure 



Figure 6. Input from Pressure Tester. 
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The response factor was computed by the following 


formula: 


AREA. 

in 


AREA 


out 

R = Response Factor 

Area. = Area under Input spectrum in stated frequency 
band 

Area^^^ = Area under response spectrum in same band 

The results of this analysis were: 

1. For frequencies greater than .15 sec"^, the response 
factor was equal to 1.00. 

2. For frequencies between .10 sec”^ and .15 sec"^, the 
response factor was equal to 1.21. 

3. For frequencies between .05 sec”^ and .10 sec”^, the 
response factor was 1.63. 

Therefore, these results suggest that the instrument over¬ 
responds to low frequency (.05-.10 sec”^) waves. Frequencies 
less than .05 sec”^ were not considered because the instru¬ 
ment wasn’t designed to measure frequencies that low. There 
is uncertainty in this analysis due to the uncertainty in 
the slope of the input pressure curve. 
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V. FIELD OPERATIONS 


This section deals with the problems and techniques 
involved in handling the instrument system at sea. 

The main difficultioss at sea occurred in handling the 
sea cable. The sea cable has a tendency to coll and kink 
whenever it is not under tension. When a severe kink or 
puncture in the insulation occurs, the sea cable does not 
conduct the signal properly. Furthermore, the sea cable 
becomes difficult to handle when wet because of its small 
diameter. In addition to these problems, the 90 lb. test 
sea cable may snap with sudden jerks. 

These problems were solved by taping a small diameter 
steel cable foot by foot to the sea cable. The rigidity 
of the steel cable prevents the sea cable from coiling and 
the tape enables the operator to get a good grip on the 
line. Furthermore, the steel cable Increases the safety 
margin against losing the transducer. 

The field operations described here were made from the 
U.S. Naval Postgraduate School’s 63 foot oceanographic 
research boat. This boat is well-suited to operations 
Involving this Instrument because the launching and recover¬ 
ing of the buoy could be done by hand. A larger vessel 
would need a water level platform or smaller boat in order 
to handle the buoy easily. It would be difficult to use 
a hoist because of the delicate nature of the instrument 
and the large antenna mounted on top of the buoy. 
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In launching the floating package, the first step is 
to secure the buoy end of the safety line. With the boat 
engines stopped, the transducer and sea cable are then 
lowered in a hand-over-hand method, care being taken not to 
allow the sea cable to scrape against anything. This is 
best accomplished with the transducer lowered on the wind¬ 
ward side of the boat. Then the sea cable and safety line 
are attached to the buoy, which is lowered by hand into 
the water. To obtain a vertical alignment with the trans¬ 
ducer, the buoy is pulled away from the boat. 

An audible tone from the receiver indicates that the 
system is in operation. It is recommended that the boat 
maintain a position 30 to 60 yards downwind from the buoy in 
order to insure strong reception of the buoy signal and to 
keep the buoy in sight. After the boat is in position and 
the equipment is functioning properly, the reference signal 
generator is turned on to indicate which portion of the 
magnetic tape has data to be processed. 

The recovery of the buoy is similar to the launching 
operation. The buoy is easily brought aboard, by hand, with 
the base of the antenna serving as a handle. The sea cable 
and transducer are then hauled slowly aboard hand-over-hand. 

These techniques for handling the instrument were satis¬ 
factory for waves as high as 12 feet with white caps present 
For higher seas, different techniques in handling the Instru 
ment may be required. 


26 












Further detail of the operation of the Instrument can 


be found In the Blssett-Berman Manual [3], 
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VI, RAW DATA 


Pour wave records were made under varying sea and wind 
conditions. All four records were made in an area centered 

10 miles WNV/ of Monterey Harbor. All measurements were 
made in water with a depth in excess of 800 fathoms. Table 

11 describes the characteristics associated with each wave 
record. 

The Douglas Sea State Code taken from Falrbrldge [4], 
describes codes 2, 3, and 4 (used in Table II) as follows: 


State 2 - small wavelets, crests of glassy appear¬ 

ance, not breaking; 

State 3 - large wavelets, crests begin to break, 

scattered whltecaps; 

State ^ - moderate waves taking longer form, many 

whitecaps, some spray. 
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DATA CHARACTERISTICS 

















VII, DATA PROCESSING 


Processing the data consisted of digitizing the wave 
record and obtaining a power spectrum. 

Because of certain characteristics In the digitizing 
process, which will be explained later in this section, 
the data on the magnetic tape were played back at a speed 
of 60 Ips (Inches per second) (recording speed was 1 7/8 
Ips) through a chart recorder moving at a speed of 
50 mm/sec. Furthermore, the data were played through a 
Khronhite filter, model 3321, with the cutoff frequency 
of 8 cps real time in order to get rid of noise introduced 
by electrical Interference from the boat. 

The wave records (Figure 8) were digitized using a 
Calma Company Model 480 Digitizer which records Ay (of 
the cartesian co-ordinates of the pressure curves) for 
every .01 Inch of-stylus travel in the x-dlrection. The 
digitized data were recorded on 5 Inch mini-reels at 556 
bytes per Inch. In order to make the data compatible to 
the IBM 360 , the record on the mini-reel was converted to 
punched cards by program "Convert" (see Appendix B) on 
the CDC 6500 at the Fleet Numerical Weather Central, 
Monterey, California. Because of a limited buffer capacity 
program "Convert" is limited to I 5 OOO sample points. There 
fore, 150 Inches is the maximum wave record length that can 
be processed. This was the reason for condensing the wave 
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Figure 8. Filtered Wave Record (major division represents two feet) 










































































































record by playing back the data at a speed increased by a 
factor of 32 times the recording speed. Under these con¬ 
ditions, one second of real data was represented by 1.56 mm 
on the chart paper. 

Due to the fixed sampling rate of the digitizer, the 
sampling interval. At, was .I 63 seconds. According to the 
Nyqulst criterion, the maximum allowable At is: 


At 


max 


1 


2f 


max 


With f equal to I .05 cps from Section IV, At = .^76 
max max 

seconds. Therefore, the actual sampling Interval, .I 63 
seconds, successfully meets the Nyqulst criterion. 

The power spectrum was obtained by utilizing IBM sub¬ 
routine RHARM. This subroutine computes one-dimensional 
Fourier coefficients by using the fast Fourier transform. 
The power spectrum was smoothed over by averaging every 20 
Fourier coefficients resulting in approximately ^0 degrees 
of freedom, The entire procedure was performed by program 
’’Spectrum” in Appendix B. 
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VIII. COMPARISON OF RESULTS WITH THEORETICAL WAVES 


Figures 9-12 in Appendix A give the power spectra of 
the four records taken by the wave meter. These spectra 
were compared with visual observaitons and with waves 
expected from storm areas indicated by weather maps obtained 
from Fleet Numerical Weather Central. 

Significant height estimates were obtained by observing 
by eye (with two different observers) and recording the 
vertical movement of the buoy from trough to crest of a 
large number of successive waves. The average of the high¬ 
est one-third of these observations gives the significant 
height. From the power spectra, the significant height was 
calculated by 

= 2.83ViT 

where 

E = 2(Variance) = area under power spectrum 

Table III and Figure 13 give the results of these computa¬ 
tions. The results suggest that the meter may have responded 
accurately to both principle spectral bands. This is not 
in agreement with the analysis concerning the known input 
and its response from Section IV. However, Figure 13 is 
only suggestive because of the. limited number of records 
taken. 

The power spectra were compared to the reports from 
the Fleet Numerical Weather Central of the combined sea 
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Number 
of V/aves 
Observed 

visual 

%/3 
(feet) 

Variance 

Meter 

«l/3 
(feet) 

May 11 

125 

5.10 

1.85 

3.85 

May 20 

75 

8.35 

9.15 

8.55 

J une 1 

288 

2.21 

l.i)5 

3 .'JO 

June 2 

26 l 

').55 

2.03 

' 1.05 


TABLE III 

RESULTS FROM SIGNIFICANT HEIGHT ANALYSIS 



Figure 13, Eye Observations vs. Meter 
Observations of Significant Height 
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height 5 which gives contour lines representing significant 
sea height. 

All the combined sea height analyses for the appropriate 
days indicate a significant height of approximately three 
feet in the vicinity of the mouth of Monterey Bay (Figures 
1^-17 in Appendix A). This roughly corresponds to all the 
reports except the one taken on the 20th of May. However, 
the combined sea height analysis may not be too reliable 
because it is based on a small number of reports for the 
area along the coast of California. 

In addition, sea surface atmospheric pressure analyses 
were used to locate and Investigate the generating areas 
of the dominant frequency bands for the swell recorded by 
the instrument. Taking the peak frequency from the power 
spectra, the group velocity of the dominant swell was com¬ 
puted using 

c = 

E 2f 

Cg = deep water group velocity in knots 
f = peak frequency in sec”^. 

Then, the distance traveled by the dominant wave was computed 
by using the difference between the time of the sea surface 
pressure analysis and the time the wave record was taken. A 
line was drawn on each sea surface pressure analysis (Figures 
18-21 in Appendix A) to represent the starting location for 
the swell observed at Monterey Bay. These lines appear to 
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lie In each case within the possible generating area of the 
dominant frequency swell. 


Computations were also made to determine whether, accord¬ 
ing to a commonly used wave generation model, winds in these 
possible generating areas produced the recorded dominant 
frequencies. The geostrophic wind speed was calculated from 
the pressure field using 



1 



AP ^ 9.33 AP 
AN sln(J) AN 


Pfipsln 


where 


V = geostrophic wind speed in knots 
S 

= angular velocity of earth 
<t> = latitude 

AP = pressure difference in millibars 

AN = latitude change (60 miles 1 degree of 


latitude) 

The surface wind speed, , is then approximated closely 
enough for present purposes by; 

V = .8V 
s g 

With values of surface wind speed and estimated fetch, the 
maximum periods which these generating areas were capable 
of producing was estimated by using co-cumulatlve spectra 
for wind speeds as a function of fetch, as given by Pierson, 
Neumann, and James [51. In all four cases, the generating 
areas apparently were capable of producing the dominant 
frequency swell recorded by the meter. Table IV gives the 
results of the computations made. 
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TABLE IV 


RESULTS FROM GENERATING AREA ANALYSIS 



Peak 
Freq. 

(sec~^) 

Period 

(sec) 

Travel 
Time 
(hrs.) 

Miles 

Traveled 

V 

s 

(knots) 

Fetch 
(miles) 

May 11 

.08 

12.5 

kl 

900 

14.7 

200 

May 20 

.11 

9.1 

35 

550 

8.2 

200 

June 1 

.08 

12.5 

65 

f. 

19*10 

9.4 

200 

June 2 

.08 

12.5 

65 

19')0 

13.0 

200 


Furthermore, the surface pressure analysis seemed to 
indicate that secondary peaks in the power spectra are due 
to local wind conditions. For example, it is noted that a 
large secondary peak not present on 1 June occurs at 
approximately .25 sec ^ on 2 June. From Section VI, it was 
recorded that there was a large Increase in the wind over 
that 2k hour period. This suggests that the secondary peak 
was due to the local winds. 


37 






IX. CONCLUSIONS 


The portable free floating wave buoy provides a rela¬ 
tively inexpensive method for gathering wave records. The 
evaluation of the records suggests that the main spectral 
bands 5 the swell and the local wind waves of the wave 
system^ are properly resolved by the meter. However^ 
further dynamic and static calibration Is suggested in 
order to resolve the discrepancy In amplitude measurements 
between the Input-response analysis and the significant 
height analysis. Comparison of the buoy’s wave records with 
the simultaneously made records of well-calibrated wave 
meters on Flip would be a desirable future project. Also, 
further comparisons with weather analyses would be profit¬ 
able In evaluation of the meter. 

The portable nature and ease of operation of the free 
floating wave buoy provide the opportunity for deep water 
wave forecast verifications to be made by organizations 
which are financially limited. The meter permits records 
of deep water waves to be taken more quickly, easily, and 
cheaply than with most other existing deep water wave 
meters. 
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Figure 9. Spectral Wave Analysis ... 11 May 70 

90^ Confidence: = 1.5S, = . 71S 
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Figure 12. Spectral Wave Analysis ... 2 June 70 

90^ Confidence: S^^ = 1.5S, S^^^ = .71S 














































Figure 1^. Combined Sea Height Analysis ... 
OOZ, 12 May 70 
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Figure 15. Combined Sea Height Analysis ... 
OOZ, 21 May 70 





































Figure 16. 


Combined Sea Height Analysis ... 
OOZ, 2 June 70 


^5 




















Figure 17 . Combined Sea Height Analysis ... 
OOZ, 3 June 70 
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Figure 18. Sea Surface Pressure Analysis 
OOZ, 10 May 70 
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Figure 19. Sea Surface Pressure Analysis ... 
06Z, 19 May 70 
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Figure 20. Sea Surface Pressure Analysis ... 
OOZ, 29 May 70 
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Figure 21, Sea Surface Pressure Analysis ... 
OOZ, 30 May 70 
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SUBROUTINE RHARM 


Given 2N real numbers, X , X, ... X„,, ,, the subroutine 

o 1 2N~1^ 

cnmputes Fourier coefficients, a2, 

^N—^N-1’ equation 


Xj - I a^ + ^a^cos(^) + bj^sin(^) 


+ I a^C-l)-’ 


where 


0 , 1 , 


2N-1 


Therefore, 
contains: 


!£ h 
n 5 t> 


with N inputs from a time series, output 


0, ^2’ ^2’ ^2^ 


^N/2’ '^N ° 


to' get the power spectrum. Program Spectrum gives the spec- 
trail estimate, S, equal to: 


|tl 1+1 r 

E S„ = E 
N=1 k=2 


+ ^k 


where T = record length. 
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